Efficient minimization of multipole electrostatic potentials in torsion space

The development of models of macromolecular electrostatics capable of delivering improved fidelity to quantum mechanical calculations is an active field of research in computational chemistry. Most molecular force field development takes place in the context of models with full Cartesian coordinate degrees of freedom. Nevertheless, a number of macromolecular modeling programs use a reduced set of conformational variables limited to rotatable bonds. Efficient algorithms for minimizing the energies of macromolecular systems with torsional degrees of freedom have been developed with the assumption that all atom-atom interaction potentials are isotropic. We describe novel modifications to address the anisotropy of higher order multipole terms while retaining the efficiency of these approaches. In addition, we present a treatment for obtaining derivatives of atom-centered tensors with respect to torsional degrees of freedom. We apply these results to enable minimization of the Amoeba multipole electrostatics potential in a system with torsional degrees of freedom, and validate the correctness of the gradients by comparison to finite difference approximations. In the interest of enabling a complete model of electrostatics with implicit treatment of solvent-mediated effects, we also derive expressions for the derivative of solvent accessible surface area with respect to torsional degrees of freedom

Protein–DNA binding specificity predictions with structural models

Protein–DNA interactions play a central role in transcriptional regulation and other biological processes. Investigating the mechanism of binding affinity and specificity in protein–DNA complexes is thus an important goal. Here we develop a simple physical energy function, which uses electrostatics, solvation, hydrogen bonds and atom-packing terms to model direct readout and sequence-specific DNA conformational energy to model indirect readout of DNA sequence by the bound protein. The predictive capability of the model is tested against another model based only on the knowledge of the consensus sequence and the number of contacts between amino acids and DNA bases. Both models are used to carry out predictions of protein–DNA binding affinities which are then compared with experimental measurements. The nearly additive nature of protein–DNA interaction energies in our model allows us to construct position-specific weight matrices by computing base pair probabilities independently for each position in the binding site. Our approach is less data intensive than knowledge-based models of protein–DNA interactions, and is not limited to any specific family of transcription factors. However, native structures of protein–DNA complexes or their close homologs are required as input to the model. Use of homology modeling can significantly increase the extent of our approach, making it a useful tool for studying regulatory pathways in many organisms and cell types

Quantitative profiling of selective Sox/POU pairing on hundreds of sequences in parallel by Coop-seq

© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research. Cooperative binding of transcription factors is known to be important in the regulation of gene expression programs conferring cellular identities. However, current methods to measure cooperativity parameters have been laborious and therefore limited to studying only a few sequence variants at a time. We developed Coop-seq (cooperativity by sequencing) that is capable of efficiently and accurately determining the cooperativity parameters for hundreds of different DNA sequences in a single experiment. We apply Coop-seq to 12 dimer pairs from the Sox and POU families of transcription factors using 324 unique sequences with changed half-site orientation, altered spacing and discrete randomization within the binding elements. The study reveals specific dimerization profiles of different Sox factors with Oct4. By contrast, Oct4 and the three neural class III POU factors Brn2, Brn4 and Oct6 assemble with Sox2 in a surprisingly indistinguishable manner. Two novel half-site configurations can support functional Sox/Oct dimerization in addition to known composite motifs. Moreover, Coop-seq uncovers a nucleotide switch within the POU half-site when spacing is altered, which is mirrored in genomic loci bound by Sox2/Oct4 complexes.Link_to_subscribed_fulltex

Rescue of Degradation-Prone Mutants of the FK506-Rapamycin Binding (FRB) Protein with Chemical Ligands

We recently reported that certain mutations in the FK506-rapamycin binding (FRB) domain disrupt its stability in vitro and in vivo (Stankunas et al. Mol. Cell , 2003 , 12 , 1615). To determine the precise residues that cause instability, we calculated the folding free energy (Δ G ) of a collection of FRB mutants by measuring their intrinsic tryptophan fluorescence during reversible chaotropic denaturation. Our results implicate the T2098L point mutation as a key determinant of instability. Further, we found that some of the mutants in this collection were destabilised by up to 6 kcal mol −1 relative to the wild type. To investigate how these mutants behave in cells, we expressed firefly luciferase fused to FRB mutants in African green monkey kidney (COS) cell lines and mouse embryonic fibroblasts (MEFs). When unstable FRB mutants were used, we found that the protein levels and the luminescence intensities were low. However, addition of a chemical ligand for FRB, rapamycin, restored luciferase activity. Interestingly, we found a roughly linear relationship between the Δ G of the FRB mutants calculated in vitro and the relative chemical rescue in cells. Because rapamycin is capable of simultaneously binding both FRB and the chaperone, FK506-binding protein (FKBP), we next examined whether FKBP might contribute to the protection of FRB mutants. Using both in vitro experiments and a cell-based model, we found that FKBP stabilizes the mutants. These findings are consistent with recent models that suggest damage to intrinsic Δ G can be corrected by pharmacological chaperones. Further, these results provide a collection of conditionally stable fusion partners for use in controlling protein stability.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/56088/1/1162_ftp.pd

Serverification of Molecular Modeling Applications: the Rosetta Online Server that Includes Everyone (ROSIE)

The Rosetta molecular modeling software package provides experimentally tested and rapidly evolving tools for the 3D structure prediction and high-resolution design of proteins, nucleic acids, and a growing number of non-natural polymers. Despite its free availability to academic users and improving documentation, use of Rosetta has largely remained confined to developers and their immediate collaborators due to the code's difficulty of use, the requirement for large computational resources, and the unavailability of servers for most of the Rosetta applications. Here, we present a unified web framework for Rosetta applications called ROSIE (Rosetta Online Server that Includes Everyone). ROSIE provides (a) a common user interface for Rosetta protocols, (b) a stable application programming interface for developers to add additional protocols, (c) a flexible back-end to allow leveraging of computer cluster resources shared by RosettaCommons member institutions, and (d) centralized administration by the RosettaCommons to ensure continuous maintenance. This paper describes the ROSIE server infrastructure, a step-by-step 'serverification' protocol for use by Rosetta developers, and the deployment of the first nine ROSIE applications by six separate developer teams: Docking, RNA de novo, ERRASER, Antibody, Sequence Tolerance, Supercharge, Beta peptide design, NCBB design, and VIP redesign. As illustrated by the number and diversity of these applications, ROSIE offers a general and speedy paradigm for serverification of Rosetta applications that incurs negligible cost to developers and lowers barriers to Rosetta use for the broader biological community. ROSIE is available at http://rosie.rosettacommons.org

Computational reprogramming of homing endonuclease specificity at multiple adjacent base pairs

Site-specific homing endonucleases are capable of inducing gene conversion via homologous recombination. Reprogramming their cleavage specificities allows the targeting of specific biological sites for gene correction or conversion. We used computational protein design to alter the cleavage specificity of I-MsoI for three contiguous base pair substitutions, resulting in an endonuclease whose activity and specificity for its new site rival that of wild-type I-MsoI for the original site. Concerted design for all simultaneous substitutions was more successful than a modular approach against individual substitutions, highlighting the importance of context-dependent redesign and optimization of protein–DNA interactions. We then used computational design based on the crystal structure of the designed complex, which revealed significant unanticipated shifts in DNA conformation, to create an endonuclease that specifically cleaves a site with four contiguous base pair substitutions. Our results demonstrate that specificity switches for multiple concerted base pair substitutions can be computationally designed, and that iteration between design and structure determination provides a route to large scale reprogramming of specificity

Redesign of the monomer-monomer interface of Cre recombinase yields an obligate heterotetrameric complex

Cre recombinase catalyzes the cleavage and religation of DNA at loxP sites. The enzyme is a homotetramer in its functional state, and the symmetry of the protein complex enforces a pseudo-palindromic symmetry upon the loxP sequence. The Cre-lox system is a powerful tool for many researchers. However, broader application of the system is limited by the fixed sequence preferences of Cre, which are determined by both the direct DNA contacts and the homotetrameric arrangement of the Cre monomers. As a first step toward achieving recombination at arbitrary asymmetric target sites, we have broken the symmetry of the Cre tetramer assembly. Using a combination of computational and rational protein design, we have engineered an alternative interface between Cre monomers that is functional yet incompatible with the wild-type interface. Wild-type and engineered interface halves can be mixed to create two distinct Cre mutants, neither of which are functional in isolation, but which can form an active heterotetramer when combined. When these distinct mutants possess different DNA specificities, control over complex assembly directly discourages recombination at unwanted half-site combinations, enhancing the specificity of asymmetric site recombination. The engineered Cre mutants exhibit this assembly pattern in a variety of contexts, including mammalian cells

Measurement of the production of a W boson in association with a charm quark in pp collisions at √s = 7 TeV with the ATLAS detector

The production of a W boson in association with a single charm quark is studied using 4.6 fb−1 of pp collision data at s√ = 7 TeV collected with the ATLAS detector at the Large Hadron Collider. In events in which a W boson decays to an electron or muon, the charm quark is tagged either by its semileptonic decay to a muon or by the presence of a charmed meson. The integrated and differential cross sections as a function of the pseudorapidity of the lepton from the W-boson decay are measured. Results are compared to the predictions of next-to-leading-order QCD calculations obtained from various parton distribution function parameterisations. The ratio of the strange-to-down sea-quark distributions is determined to be 0.96+0.26−0.30 at Q 2 = 1.9 GeV2, which supports the hypothesis of an SU(3)-symmetric composition of the light-quark sea. Additionally, the cross-section ratio σ(W + +c¯¯)/σ(W − + c) is compared to the predictions obtained using parton distribution function parameterisations with different assumptions about the s−s¯¯¯ quark asymmetry

Measurement of the angular coefficients in Z-boson events using electron and muon pairs from data taken at √s=8 TeV with the ATLAS detector

The angular distributions of Drell-Yan charged lepton pairs in the vicinity of the Z-boson mass peak probe the underlying QCD dynamics of Z-boson production. This paper presents a measurement of the complete set of angular coefficients A0−7 describing these distributions in the Z-boson Collins-Soper frame. The data analysed correspond to 20.3 fb−1 of pp collisions at s√=8s=8 TeV, collected by the ATLAS detector at the CERN LHC. The measurements are compared to the most precise fixed-order calculations currently available (O(α2s))(O(αs2)) and with theoretical predictions embedded in Monte Carlo generators. The measurements are precise enough to probe QCD corrections beyond the formal accuracy of these calculations and to provide discrimination between different parton-shower models. A significant deviation from the (O(α2s))(O(αs2)) predictions is observed for A0 − A2. Evidence is found for non-zero A5,6,7, consistent with expectations

Measurement of the low-mass Drell-Yan differential cross section at √s = 7 TeV using the ATLAS detector

The differential cross section for the process Z/γ ∗ → ℓℓ (ℓ = e, μ) as a function of dilepton invariant mass is measured in pp collisions at s√ = 7 TeV at the LHC using the ATLAS detector. The measurement is performed in the e and μ channels for invariant masses between 26 GeV and 66 GeV using an integrated luminosity of 1.6 fb−1 collected in 2011 and these measurements are combined. The analysis is extended to invariant masses as low as 12 GeV in the muon channel using 35 pb−1 of data collected in 2010. The cross sections are determined within fiducial acceptance regions and corrections to extrapolate the measurements to the full kinematic range are provided. Next-to-next-to-leading-order QCD predictions provide a significantly better description of the results than next-to-leading-order QCD calculations, unless the latter are matched to a parton shower calculation